Fluctuation theorem for currents and nonlinear response coefficients 
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We use a recently proved fluctuation theorem for the currents to develop the response theory 
of nonequilibrium phenomena. In this framework, expressions for the response coefficients of the 
currents at arbitrary orders in the thermodynamic forces or affinities are obtained in terms of 
the fluctuations of the cumulative currents and remarkable relations are obtained which are the 
consequences of microreversibility beyond Onsager reciprocity relations. 



I. INTRODUCTION 
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Onsager's classic work of 1931 [1J has shown that the linear response coefficients relating the currents to the ther- 
modynamic forces or affinities obey reciprocity relations as a consequence of the reversibility of the underlying 
microscopic dynamics. Another consequence of microreversibility are the so-called fluctuation theorems, which charac- 
terize the large fluctuations of physical quantities in nonequilibrium systems. They have been derived in deterministic 
IHIEU or Markovian stochastic systems 0, H H, E3] 

and concern different quantities such as the entro py p roduction 
^ . I], the dissipated work [HI, [Hj], or the currents crossing the system in a nonequilibrium situation (l3.1l4l [l5|. Such 

S relations are important because they are valid far from equilibrium. Close enough to equilibrium where the response 
of the system is linear in the affinities, the Onsager reciprocity relations can be deduced from the fluctuation theorem 



On the other hand, it is known that far-from-equilibrium systems may present nonlinear responses to nonequilibrium 
^ ■ constraints. The response is said to be nonlinear if the currents crossing the nonequilibrium system depend nonlinearly 
on the affinities. The coefficients characterizing such nonlinear responses can be obtained by expanding the currents 

S in the powers of the affinities. The terms linear in the affinities are the linear response coefficients obeying Onsager's 
reciprocity relations. The terms which are quadratic, cubic, quartic, etc... in the affinities are called the nonlinear 
, response coefficients. We may wonder if the nonlinear response coefficients would obey relations beyond Onsager's 
£^ ■ ones as the conequence of the fundamental microreversibility. 

The purpose of the present paper is to show that, indeed, the nonlinear response coefficients do obey remarkable 
relations which have their origin in microreversibility. For this purpose, we use a fluctuation theorem for the currents 
which was first proven for mechanically driven Markovian processes |8|, I fA 1 1011 , then in the more general framework of 
7— I ■ Schnakenberg network theory [l6[ which includes reactive processes jl3l. Il4f . as well as in non-Markovian situations 
[l7| . This fluctuation theorem directly concerns the generating function of the different fluctuating currents crossing a 
nonequilibrium system. Consequently, the nonlinear response coefficients can be directly obtained from the generating 
function by successive differentiations, so that the symmetry of the fluctuation theorem for the currents can be used 
in a straightforward way. The fluctuation theorem for the currents has been proved elsewhere [14j under the general 
conditions enunciated by Schnakenberg [l6| and we start from this important result to obtain remarkable relations as 
| the consequences of microreversibility. 

The plan of the paper is the following. In Section |TT1 we summarize the results about the fluctuation theorem 
for the currents. Section IIIII is devoted to the derivation of the consequences of the fluctuation theorem on the 
response coefficients up to the cubic response coefficients with comparisons with known results. In Section IIV1 the 
nonlinear response coefficients are systematically calculated and the generalizations of Onsager relations are obtained 
at arbitrarily large orders. Conclusions are drawn in Section IVl 



II. FLUCTUATION THEOREM FOR THE CURRENTS 

The fluctuation theorem for the currents relates the probability of observing given values for the cumulated currents 
to the probability of observing negative values via the following exponential relation valid in the long-time limit: 
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where j 7 (£) denote the independent fluctuating currents and A~, are the corresponding affinities (also called the 
thermodynamic forces) driving the system out of equilibrium 0, [l6[ • 
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If we introduce the decay rate of the probability that the cumulated currents take given values 



1 



H({^}) = lim InProb 



t-»oo t 

the fluctuation theorem can be written as 

h({-^})-h({^}) = J2M 



dt'j^t') g (£ 7 ,£ 7 +d£ 7 ) 



(2) 



(3) 



The Legendre transform of the decay rate ([2|) 



Q({A 7 }) = Max {M 



is the generating function of the currents defined by 



#({£ 7 }) + E A ^ 



1 



Q({A 7 };{A 7 }) = lim --ln(e 



t^oo t 



E T A 7 G 7 (t) 



in terms of the cumulative currents also called the Helfand moments [181 ] : 

G 7 (t)= [ j-y(t')dt' 

Jo 

The fluctuation theorem for the currents ([T]) is now expressed as 

Q({A 7 };{A y }) = Q({A y -A 7 };{A y }) 



(4) 



(5) 



(6) 



(7) 



in terms of the generating function. This relation has been derived in the context of stochastic processes [TJ]- In 
this description the system is described by a probability distribution over the possible states and which is ruled by a 
master equation with several transition rates. The macroscopic affinities A 1 are then identified using Schnakenberg's 
network theory (l6j . 

In the nonequilibrium steady state, the mean value of the current j a (t) is given by 



Jot 



dQ 



= lim - / {j a (t'))dt' = lim ~(G a (t)) 

{A 7 =0} t^oo t Jq t^oo t 



(8) 



III. CONSEQUENCES FOR NONLINEAR RESPONSE 

In this section, we prove that the fluctuation theorem ([7| for the macroscopic currents ((8|) have important conse- 
quences not only on the linear response coefficients but also at the level of the nonlinear response. In general, the 
macroscopic currents can be expanded as power series of the macroscopic affinities: 



Ja = ^2 LapAp + - ^2 MafiyApAj + - ^ N^gApA^As 



(9) 



/3,7 



0,7,6 



The linear response of the currents J a with respect to a small perturbation in the affinities Ap is characterized by the 
Onsager coefficients L a p, and the nonlinear response by the higher-order coefficients M a p-y, N a p-ys,... 



A. Onsager reciprocity relations 

The Onsager coefficients are defined close to the equilibrium in terms of the generating function J7J) by 

d 2 Q 



dJ a 



J af3 



dA t 



A=o d\ a dA, 



-(0;0) 



(10) 
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If we differentiate the expression ([7]) of the fluctuation theorem with respect to X a and Ap we find that 

d 2 d 2 d 2 

' (A;A)=-^-(A-A;A)-^ r -(A-A;A) (If) 



8X a dAp ' d\ a d\p ' d\ a dAp 

Setting A = and A = 0, we obtain the relation 



2 O d 2 

" -(0;0) = -^f^(o ; o) (12) 



dX a dAp ' d\ a d\fj 
or 

^--bSb** ' (13) 

as already shown in reference [U, Q . Hence the Onsager reciprocity relations 

L a /3 = Lp a (14) 

We notice that no further relation is obtained by differentiating the fluctuation relation twice with respect to 
either the parameters A or the affinities A. 

B. Green-Kubo and Einstein-Helfand formulas 

By using Eq. I|13p. we obtain the Onsager coefficients as 

Lap = \ f +OC ([j a (t) - (j a )] fo(0) ~ 0»]) eq dt = lim 1 (AG Q (t)AGp (t)) cq (15) 

in terms of the time correlation functions of the instantaneous currents or the corresponding Helfand moments: 

AG a (t) = G a (t)-(G a (t)) (16) 

Here, the statistical average is carried out with respect to the state of thermodynamic equilibrium. In Eq. (|15p . 
the formulas giving the coefficients in terms of the time correlation functions are known as the Green-Kubo formulas 
[l^.[20| (or the Yamamoto-Zwanzig formulas in the context of chemical reactions [2l|,[22j]). The other formulas giving 
the coefficients in terms of the Helfand moments or cumulative currents are known as the Einstein-Helfand formulas 
11. 



C. Relations for the second-order response coefficients 

The second-order response coefficients are defined as 

M ^ E ^ (0;0) (17) 

in terms of one derivative with respect to the parameter X a generating the current J a and two derivatives with respect 
to the affinities Ap and A 1 . 

Our purpose is to relate these nonlinear response coefficients to quantities with a reduced number of derivatives 
with respect to the affinities, thus characterizing the fluctuations instead of the response. 

Such relations are obtained by continuing the procedure started to get the Onsager reciprocity relations by further 
differentiating the generating function. If we differentiate the identity (jlip with respect to A 1 and set A = and 
A = 0, we obtain the second-order response coefficients as 

Ma ^ = -LaJa^a 7 (0;0) - LaJa^a 7 (0;0) - LaJva/ 0;0) (18) 

Using the symmetry ([7]) of the current fluctuation theorem at equilibrium, we see that the generating function 
Q(A; 0) = Q(— A; 0) is an even function of A 7 at equilibrium. The first term of the right-hand side of equation (fT8|) is 
a third derivative with respect to the parameters A 7 calculated at equilibrium, hence it must vanish identically. 
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by 



The other terms in Eq. (|18p are related to the spectral function of the nonequilibrium current fluctuations defined 



:[j«(t) - Ua)} fo(o) - o»]> dt 



(19) 



where the statistical average is here taken with respect to the nonequilibrium steady state. Here, we introduce the 
quantities 



Ra/3,' 



d 3 Q 



dXadXpdA-f 
q r+oo 



(0;0) 



dA 



([?*(*) ~ (j a )} b>(0) - (?»]> ^ 



7 J — oo 

' } lim \(AG a (t)AG p (t)) 



A=0 



dAj t^oo t 



A=0 




dA, 



■S a/3 (w = 0) 



A=0 



(20) 



which characterize the sensitivity of the current fluctuations out of equilibrium. Equation (|20p shows that the sensitiv- 
ity coefficients are given in terms of the derivative with respect to the affinities of the spectral function or, equivalently, 
of the diffusivities of the nonequilibrium currents defined by 



A=0 



(21) 



According to the fluctuation theorem and Eq. (fT5)) . we find that the second-order response coefficients are given in 
terms of the sensitivity coefficients by 



For the case j3 — 7 we find 



M, 



a/3/? 



R 



a/3,/3 



(22) 



(23) 



In particular the response coefficients M a ^ present the expected symmetry M a py = M a ~ f p. The second-order 
coefficients are thus related to the diffusivities by 



M, 



oi/3-y 



d 



dA, 



D 



8 



a/3 



dAp a \ 



(24) 



A=Q 



Thanks to the fluctuation theorem for the currents, we can therefore relate the second-order nonlinear response coef- 
ficients to quantities characterizing the nonequilibrium fluctuations such as the spectral functions or the diffusivities 
of the currents in the nonequilibrium steady state. We notice that the number of derivatives with respect to the 
affinities has indeed been reduced. 

Similar expressions can be found for even higher-order relations where the odd derivatives with respect to the A's 
automatically vanish at equilibrium. 



D. Relations for the third-order response coefficients 



A similar reasoning can be carried out for the third-order response coefficients defined by 

d 4 Q 



N. 



-(0;0) 



dXadAfidAydAs 

Differentiating the identity (|11|) twice with respect to the affinities A 1 and As shows that 
1 



Naff^S = — 



d 4 Q 1 1 

(0; 0) — - (S a p lt s + S a ps,i + S ai sj3) + 77 (TaPrfS + T ajt ps + T a s,i3i) 



2 9A Q 9A (3 9A 7 9A < 5 



(25) 



(26) 
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with 



and 



T, 



a/3, 7(5 



d 4 Q 
dX a dXpdX 1 dAs 



d 4 Q 



(0;0) 



-(0;0) 



dXadXpdAjdAg 

Differentiating Eq. (jTTJ) twice with respect to the parameters A 7 and Xs shows that 

f d 4 Q 



S a f3-y,8 



-(0;0) 



2 dXadXpdX-fdXs 

which proves the total symmetry of this tensor. Accordingly the third-order response coefficients are given by 



We thus obtain the reciprocity relations that the fourth-order tensor 

2N a j3jS — T a p^s — T al .ps — T a s t f3 7 

must be totally symmetric. 

The tensor (|27p can be expressed as 



^- lim ^(AG a (t)AG /J (t)AG 7 (i)> 



,4=0 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 



which do not vanish in general as for the even-order cases. Equation (|3"2"|) shows that the tensor (|2T|) characterizes 
the sensitivity of the third-order moments of the cumulative currents with respect to the noncquilibrium constraints. 
Moreover, the expression ([29]) shows that this tensor can also be calculated at equilibrium as 



S a /3j,S 



lim 

t — >oo 



l r 

2t 



(AG a {t)AGp (t) AG 7 (t)AGj (t)) 

(AG a (t) AG/3 (t) ) (AGy (t) AGs (t)> 
(AG a (t) AG 7 (t) } (AG (t) AGs (t)> 

( AGq, (t) AG« (t) } ( AG^g (i) AG 7 (i) ) 



(33) 



where it characterizes the fluctuations. The equality between Eqs. 
of the fluctuation theorem. 

On the other hand, the tensor (|28p is given by 



and 



is another remarkable consequence 
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d 


dA 1 


OAs 


d 


d 


dA 7 


dA 5 


d 


d 


dA 1 


OAs 



( \ja(t) ~ (J'a)] [7>(0) - (?»]> dt 



A=0 



S Q( 3(w = 0) 



A=0 



(34) 



Accordingly, the tensor (|28p also characterizes the sensitivity of the nonequilibrium fluctuations but now in terms of 
the second derivatives of the power spectrum with respect to the affinities. We notice the similarity with Eq. (I20[) . 
Again, the number of derivatives with respect to the affinities has been reduced compared to the definition (f25l) of 
the third-order response coefficients and this thanks to the fluctuation theorem for the currents. 
The expansion can be carried out to higher orders as done in the next section. 
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IV. RELATIONS AT ARBITRARY ORDERS 



In a macroscopic description, we consider general affinities A a conjugated to currents J a . The mean value of the 
currents can be developed as 



E^i E ^l.A,...A (35) 



=1 " A-,,. 



with the coefficients 



C 



d n J a 



(n) = 

al3 -^ dA ---dA„ 



d n+1 Q 



A=a dXadA/3 ■ ■ ■ dA^ 



(36) 

X=0,A=0 



contain n+l indices. The expansion in powers of the affinities gives Onsager's linear response coefficients as well 
as higher-order coefficients C^ n ' characterizing the nonlinear response of the system with respect to the nonequilibrium 
constraints {A € }. 

We now want to use the fluctuation theorem for the currents J7J to obtain expressions for the response coefficients. 
To do so we will consider expressions for the derivatives of Q at arbitrary orders. This will provide us with several 
non trivial relationships and we will have to combine them to obtain a simple form for the response coefficients. 

Using the fluctuation theorem (J7]), the derivatives of Q are given by 

n—k 

q^„.. ct = (~i) k £ C5;L }p ( 37 ) 

where the derivatives are calculated at A = A = 0. The notation Qa-'»rj,p---o- means that we have taken n derivatives 
with k of them corresponding to Aa,...^ and n — k corresponding to A p ,...,A a : 

= d\ a --jdt dA p .-^M~^ 0) (38) 

k n — k 

These derivatives correspond to the response of the k th cumulant with respect to the macroscopic affinities: 

Qi k :t P -« = (~l) fc+1 lim \{{G a ---G v )) (39) 

oA p ■ ■ ■ oA a t-*oo t A=0 

where ((•)) denotes the cumulant. The notation {-} p means the symmetrized ensemble with the derivatives taken 
after the term p. For example, {a/3jS}i = a, /3~/5 + j3, ajd + 7, a/36 + 6, afij. In the same way, {a/3^6}2 = a/3, j6 + 
cvy, (36 + a6, 7/? + (3"f, a6 + (36, aj + j6, a(3. There are m\/[p\(m — p)l] numbers of terms if there are m terms in the 
ensemble. The derivative (|3T)l is thus expressed as the sum of 2 n ~ k terms. The sign (— l) fc comes from the derivatives 
with respect to A while the structure of the sum comes from the derivatives with respect to A. Indeed, each such 
derivative generates two terms, one with a derivative with respect to A and the other with respect to A as can be seen 
from Eq. (0. Noting that the term p = in Eq. (|3"T)) is the same as the left-hand side of the equation, we have 

n — k 

[1 + (-i) fe+1 ] c£''t-« = (-i) fe E Q^nit^ ( 4 °) 

p=l 



so that 



= E Q^nit^ k even ( 41 ) 
P =i 

7i — k 

2 Qt'to-* = E QZt%L }p k odd (42) 

p=l 

As explained above each derivative Q^ k ' n ^ can be expressed as derivatives with respect to n — k affinities of a cumulant 
of order k at equilibrium. We thus have a number of non trivial relations between different moments and their 
derivatives with respect to the affinities calculated at equilibrium. 
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In particular, the response coefficients (|36|) of order n—1 are given by 

n— 1 

f,{n—X) _ n (l,n) _ ST^ n( 1+p '™) (A1\ 

^aP-o — ^a,/3-a- ~ 2 ^ V a{/3---er} p 

where we used Eq. (|42|) with fc = 1 which is odd. The response coefficients are thus expressed as a sum of 2 ,l_1 — I 
terms with n — 1 different tensors. 

The relations (|37|) can be used to simplify the expressions of the response coefficients. For example, the relations 
(|II j) -([4"2 f with k = n - 1 give 

n (n,n) _ 



n odd (44) 

I 



2 QL"^ = - Q^-t n even (45) 



so that the derivatives with respect to A vanishes if n is odd and we obtain the total symmetry of the tensor Q^.-r/'P 
if n is even: 



d ,. I 



1 .. I 



lim -«G tt ---G„)) = - lim -UG a ---G n G CT ))^ (46) 



which is totally symmetric. These are non-trivial consequences derived from the fluctuation theorem. If k = n we 
do not get any new informations. If k — we find a constraint on the sum of the response coefficients but it can be 
recovered from their expressions (|43[) . In fact all relations (|4f [) with fc even can be recovered from relations (|42[) for k 
odd. We now have to use them to simplify the expression of the response coefficients. 

To do this, we will thus choose to use relations (|42| to eliminate all terms of the form Q' fe > n ) where k is odd. 
However, the terms in the right hand side of Eq. (|43p are all totally symmetrized wrt the n — 1 indexes {(3 ■ ■ ■ a} while 
relations (|42p are not. The first step is thus to symmetrize relations (j4"2"|) to get 

- 4 £ ( * - 1 ) «8£s>. for fc ° dd (4?) 



Z=fc 



where 



I \ U 



(48) 



These coefficients are obtained by symmetrizing the relations (|4"!2)) . As all relations become thus totally symmetric, 
there are ^ ^ ' ^ ^ identically terms Q^g.'.^j., arising by this procedure. These coefficients are given by the numbers 

of terms in Q^,' 7 ^ , times the number of terms in Q^ +1 ' r ™- ) , divided by the number of terms in Q^ta i ■ 

Unfortunately expression (|4T|) is expressed in terms of Q^ k,n ^ with k odd and even. We will thus have to use this 
relation recursively to eliminate all terms with k odd within itself. 

One have then to eliminate successively the terms Q, 7 ^ x , Q ,a i ,—,Q^fa i 

When this is done, one can express the response coefficients in the form 

c&2 = 7i ogSU + 73 g«$L }l + • • • + 7n- 2 < n odd ) ^ 49 ) 

and 

ci;: 1 ] = 7i qS?.. ct}i + 73 q ( :i;L }3 + ■ ■ ■ + Tn-x (» even ) ( 5 °) 

where we used that Q( n > n ) = for n odd. We find here the important property that the coefficients ( ^ ' ^ ^ are 

independent of n which implies that the coefficients 7, do not depend on n. The odd response coefficients are thus 
expressed in terms of (n — 1)/2 different tensors which is far better than the n—1 tensor needed in expression (|43[) . The 
even response coefficients are thus expressed in terms of n/2 different tensors. By construction, they are symmetric 
for the permutations of the n — 1 indexes {/?••• a} as it should. 
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In particular, the tensor 



71 Q« {/S :.. ff}l + 73 Q Q{/5 ... CT}3 + • • • + 7n-3 Q{/3 ... CT}r , 



7n-l QiTt (51) 



with n even is totally symmetric. 

We now want to calculate the coefficients 7.;, with i odd, associated with the terms Q^fp^x.- The first term 71 

takes the value —1/2. The next ones are given by the successive elimination of the terms Q^fp,..^^ Q^a{p---a}i^'^ 

Q^/s 1 . Each successive elimination will change the coefficients in front of the Q ta x ■ We thus introduce 

numbers \k which denote the coefficients pondering the terms Q^^g.'.?/)-^ at the fc th successive elimination. From Eq. 

»{/3--ct} 2 



43|) we set xd = — 1/2 V/. The first step is thus to eliminate the term Q^?g... a \ so that 




*=< v ; i=2 (52) 



X& iH = l 

according to Eq. (|47j) . We can continue and eliminate the term Q a !p... lT x 4 to get 




and after fc steps we have 



Xi-|l ; )xt if^>5 
X 2 =<! V / (53) 

M if z = 4 

, xi otherwise 



Xl-H ^ (54) 




.xLi 



otherwise 



These numbers are independent of n as it should. The coefficients 7$ in Eq. I|67p - (|68p are then given by x\ once they 
remain invariant that is when A: > (i — l)/2. This construction can be summarized in the form : 

*-jE&(,li)-5 < 55 » 

where the notation ^2p= a means that we sum from p = a to b by step of 2. We also absorbed in the expression of 
7, the factor — ^ in front of expression (|43]) . The constant —\ is the contribution from the original terms Q( I+1 ' n ) in 
Eq. (|43[1 . The coefficients £ p give the number of terms Q^{j3... a } ! coming from the elimination of the previous odd 
terms. They are given by 

^=1-5 e fen (56) 

(=3 

We then find 



£3 = 1, £5 = -2, £7 = y , £9 = -62, . . . (57) 



Injecting those numbers in relations (|55|) we find 



11 1 17 31 

7i = 73 = 4, 75 = -y 77 = y, 79 = -y, ■ ■ ■ (58) 
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A shorter relation in order to obtain the ji , i > 1, is given by 



1 1 i_2 ' 2 

= --- F 

4 9 ^ 



Z=3 



i - 1 
/ - 1 



(59) 



We have thus obtained expressions for the response coefficients of order n in terms of microscopic correlation 
functions and their response to the affinity. Using the symmetry of the current fluctuation theorem we were able to 
simplify the original expressions (|4"5)) . We found a simple structure in terms of (n— l)/2 or (n/2) different tensors if n 
is respectively odd or even. These tensors characterize the fluctuations of the currents and their response to affinities 
up to an order inferior to the response coefficient. The numerical coefficients pondering the different tensors turned 
out to be independent of n, so that the same expressions arise independently of n and of its parity. Nevertheless, a 
difference arises between odd and even response coefficients: for n even there exists a totally symmetric part arising 
in the expression of that does not appear for n odd. 

For example, using (j6"8j) Onsager's coefficients are given by 



Ci'i = 7i 



)(2,2) 



)(2,2) 

la/3, 



= L 



a/3 



(60) 



according to Eq. (fTS"]) in terms of the time correlation functions of the instantaneous currents [n| [2(| or the corre- 
sponding Helfand moments [18|. Onsager's symmetry [l| is therefore verified. Here, the statistical average is carried 
out with respect to the state of thermodynamic equilibrium. 
The second-order response is expressed as 



a 



(2) 
a{3y 



n( 2 ' 3 ) 

71 Va^h 



f 



,(2,3) 



2 ^a{/3 7 }i 



}(2,3) 
" a/3,7 



which gives the response coefficients in terms of the expressions ([20| : 



i(2,3) 
■* a/3,7 



-R, 



a/3,7 



The third-order response coefficients are given by 



,(2,3) 



(61) 



(62) 



r (3) _ n (2,4) n (4,4) 
U a/?7S — 71 t ^ Q { / 3 7 «5}i + 73 ^ Q/ 3 7 5 



f r 

2 



o (2 - 4) 



^ a 



4,4) 
73 7 <5 



(63) 



where Q^ 4,4 ) is the cumulant of order 4 calculated at equilibrium by Eqs. (f3"3")) and Q*- 2 ' 4 - 1 are the second derivatives 
of the power spectra with respect to the affinities at equilibrium by Eq. (134)) : 



(2A) 



2<S'a/3 7 5 
"2~ck/3, 7 5 ^-ck 7 ,/35 ^a<5,/3 7 



(64) 
(65) 



To illustrate the elimination of the Q^ k ' n ^ with fc even in a non trivial case, let us consider the fourth-order response 
coefficients. They are given by 



C. 



(4) 

a./3-ySe 



i 4 



(1+P,5) 
a{/3 7 5e} p 



(2,5) 

a{/3 7 (5e}i 



o (2>5) 



£,(3,5) 

1 o (4 ' 5) 

2 ^ Q {/3 7 <5e} 3 



j(4,S) 
'a{/3 7 i5e}3 



(66) 



where we used that Q a p-ySe. = to go from the first line to the second and the relations (|4"2"|) to go from the second to 
the third line. 

The fourth-order response coefficients can thus be expressed in terms of the third derivatives of the spectrum with 
respect to the affinities and in terms of the first derivative of the fourth-order correlation functions with respect to 



the affinities. In particular, it presents the expected symmetry C 



(4) 
/3 7 <5e 



c 



(4) 



= •■■ = c 



(4) 



The 
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expression (|66|) could have been obtained immediately using the general form (|67| with coefficients % given by (j58[) 
One can also note that the coefficients 7$ are the same as in expression (|5T|) as it should. 

In the same way, the fifth-order response coefficients is immediately given by Eq. (|68p which reads 



C 



(5) 

a./3'ySefi 



1 r 

2 



)(2,6) 

'a{^7(5e/i}i 



)(4,6) 

'a{/375e^} 3 



)(6,6) 



as can be verified using using relations (|45f and (|42]) on the expression (|43|) of the tensor. 



Eventually, we can construct in the same way the higher-order relations for the fluctuations. Indeed the reasoning 
remain unchanged when considering the fluctuations and their responses. Using relations (|42p recursively yields, for 
odd m, 

✓~«(ro,n) ^-.(m+l.n) ^(m+3,n) ^{n—l.n) , / an\ 

Q ai - am J3-a = 71 Q ai ... a J{p...a} 1 + 73 Q ai - a J{0-a} s +"' + 7n-2 O ar .. am{(3 ... ff }„_ m _ 2 (n odd) (67) 

and 

Q ai -L m ,p- a = 7i Q„ 1 ... a J{/j...a} 1 + 73 Q ( ai ... a J {p ... a}3 + ■ ■ ■ + ln-x Q ai ... am p-a (« even) (68) 

The fluctuations are thus expressed in terms of the independent tensors with the same ponderation as for the response 
coefficients. 



V. CONCLUSIONS 



In this paper, we have shown that the fluctuation theorem for the currents (JTJ) or Q implies not only Onsager's 
reciprocity relations [l[ along with the Green-Kubo and Einstein-Helfand formulas H, [l^, [2(| [23| for the linear 
response coefficients, but also further remarkable relations for the nonlinear response coefficients at arbitrarily high 
orders. These results find their origin in the validity of the fluctuation theorem for the currents far from equilibrium 
in stochastic rate processes. The obtained relations are thus the consequences of the microreversibility. 

The response coefficients are defined by expanding the currents crossing the nonequilibrium system in powers of the 
affinities (or thermodynamic forces). Therefore, the response coefficients are defined with respect to the equilibrium 
state where the affinities vanish. Nevertheless, we can estimate the currents further away from equilibrium if we 
use an expansion up to high powers of the affinities. This explains that we need a general property valid far from 
equilibrium, such as the fluctuation theorem for the currents, in order to obtain relations on the nonlinear response 
coefficients at arbitrary orders. 

Typically, the relations described in the present paper connect quantities measuring the statistical correlations 
among m fluctuating cumulative currents to corresponding quantities among m — 1 of them with an extra derivative 
with respect to an affinity. The former characterizes the fluctuations at 771 th order and the latter the sensitivity of 
the fluctuations at the lower (777 — l) th order with respect to the nonequilibrium constraints. This is the case for 
instance for the equality between Eq. (|32[) which measures the sensitivity of the nonequilibrium correlations among 
three cumulative currents under changes of an affinity and Eq. (|33p which characterizes the fluctuations between four 
cumulative currents in the equilibrium thermodynamic state. This is the feature which is found in the relations here 
described. 

In conclusion, the theory developed in the present paper provides a general framework to formulate the nonlinear 
response theory in nonequilibrium processes. The results here reported have important applications for the nonlinear 
response properties of many nonequilibrium systems such as the chemical and biochemical reactions [l3j . the full 
counting statistics in mesoscopic conductors 15], the effusion of ideal gases [24]], and Van den Broeck's demons [25| . 
The present results could be especially important in nonequilibrium systems at the micro- and nano-scales, where the 
nonlinear response properties turn out to be dominant [261 ] . 
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